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We study the effect of thermal and quantum fluctuations on the dynamical response of a one¬ 
dimensional strongly-interacting Bose gas in a tight atomic waveguide. We combine the Luttinger 
liquid theory at arbitrary interactions and the exact Bose-Fermi mapping in the Tonks-Girardeau- 
impenetrable-boson limit to obtain the dynamic structure factor of the strongly-interacting fluid at 
finite temperature. Then, we determine the drag force felt by a potential barrier moving along the 
fluid in the experimentally realistic situation of finite barrier width and temperature. 


I. INTRODUCTION 

Superfluidity is one of the most dramatic manifes¬ 
tations of quantum many-body physics. The question 
whether a degenerate neutral bosonic fluid should be con¬ 
sidered as superfluid or not is far from trivial, especially 
in one dimension, and it is agreed that superfluidity is 
rather a collection of phenomena than a well-defined phe¬ 
nomenon. In the last decades, one-dimensional (ID) sys¬ 
tems, in which the effect of interactions is known to be 
enhanced compared to 3D, have attracted an increasing 
interest In particular, cold atomic gases can now 

be confined to investigate the effects of low dimension¬ 
ality OIMD]- An issue is the effect of the interactions 
in the degenerate fluid on the possibility of a superfluid 
flow. While several criteria lead to the same conclu¬ 
sions in 3D EHm, various aspects of coherence and su¬ 
perfluidity may be characterized in one-dimensional sys¬ 
tems by different observables. For instance, there is no 
Bose-Einstein condensation in ID owing to the Mermin- 
Wagner theorem, yet there may still subsist off-diagonal 
quasi-long-range order naiin]. 

Among other possibilities, superfluidity can be charac¬ 
terized by the absence of a (viscous) drag force acting on 
a moving fluid when it encounters an impurity or a poten¬ 
tial barrier. The case of a mobile impurity has been the 
object of intense activity in the last years [2QH26] . Here, 
we focus on an external potential barrier driven at con¬ 
stant velocity across the fluid. In this configuration, the 
drag force concept has already been used to probe super¬ 
fluidity in a 2D Bose gas m With an idealized model 
of a delta-potential, Pitaevskii and Astrakharchik [28] 
showed that, according to the drag force criterion, quan¬ 
tum fluctuations give rise to a breakdown of superfluidity 
at large interactions in a ID Bose gas, while the fluid may 
well exhibit a behavior close to superfluid when interac¬ 
tions are small. The calculation of the drag force in linear 
response theory requires the knowledge of the dynamic 
structure factor. This itself is an important measurable 
quantity which shows the many-body spectrum of collec¬ 
tive excitations in the fluid, and is experimentally acces¬ 
sible by Bragg scattering P [29ll34] . 

In this work we study the dynamic structure factor and 
the drag force under the experimentally relevant condi¬ 


tions of finite temperature and finite width of the barrier. 

We describe the system by the Lieb-Liniger model of 
interacting ID bosons [35], now realizable with exper¬ 
iments with ultracold atomic gases in a wide range of 
interaction strength [SSI EH. We focus on the strongly 
interacting regime and describe it using two complemen¬ 
tary techniques: the Tomonaga-Luttinger liquid (LL) 
and the exact Tonks-Girardeau solution (TG). The LL 
approach is valid at low energies and temperatures, for 
intermediate to strong interactions. The TG exact solu¬ 
tion describes impenetrable bosons with infinite repulsive 
interactions by means of a Bose-Fermi mapping. 

First, we determine the dynamic structure factor us¬ 
ing the LL and TG methods, thus obtaining its tem¬ 
perature dependence for an interaction regime where the 
Bogoliubov approximation [38| is not applicable. Our ap¬ 
proach is particularly suitable to describe the umklapp 
region of momenta q in the vicinity of 2k where kp 
is the Fermi wavevector of the mapped Fermi gas. In 
this region, a Random Phase Approximation perturba¬ 
tive approach [39ti42] is not applicable and predicts a 
non regular behavior for the dynamic structure factor. 
Although the dynamic structure factor is amenable to 
Bethe Ansatz calculations gSlIll], to date the strongly 
interacting regime has not been covered by this technique 
at finite temperature. In our work, specifically we focus 
on the comparison between the results obtained with the 
LL and TG approaches. In order to achieve this goal 
we determine the parameters of the LL model by solving 
the Bethe Ansatz equations for the Lieb-Liniger model 
both at zero and finite temperature. Then, in a linear 
response approach, we determine the drag force behav¬ 
ior as a function of the barrier velocity, both at zero and 
finite temperature, including the effect of a finite width 
of the barrier. 

The paper is organized as follows: in Section [IT] we 
present the system and the concepts of dynamic structure 
factor and drag force; in Sections [TTI| and [TV| we evaluate 
them for a Tonks-Girardeau gas and a Luttinger liquid 
respectively. 
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II. THE SYSTEM, AND CONCEPTS TO PROBE 
ITS SUPERFLUIDITY 


We consider a system of Nq ultracold spinless interact¬ 
ing bosonic atoms confined in a one-dimensional waveg¬ 
uide, as described by the Hamiltonian 


Ho = 



—f dx'V(x — x')n{x)n{x' 
2m 2 Jq 


( 1 ) 


where L is the length of the system, m the mass of an 
atom, -0 a bosonic field operator satisfying the commu¬ 
tation relation ['^(x),= 8{x — x'), and n{x) = 
'0'*' (x)'0(x) the density operator. We describe interactions 
as a zero-range pairwise interaction potential V{x — x') = 
g6{x — x') [45] as in the Lieb-Liniger model, with ^ > 0 to 
describe repulsive interactions. We will assume that the 
system is homogeneous in the longitudinal direction and 
adopt periodic boundary conditions. We focus on the 
concept of superfluidity as probed by a weak potential 
barrier stirred in the Bose fluid. A possible experimen¬ 
tal setup for this barrier is a Gaussian laser beam whose 
waist w moves at a velocity v in the frame of the fluid, 
leading to a perturbing Hamiltonian 


I 2 Uh 2{x — vt)‘^ , 

Hi = dx\ - -e ' 0 '*'(x)' 0 (x), ( 2 ) 

Jo V TT ^ 

Factors have been chosen so that, in the limit re ^ 0, 
we recover the results found in [28| for a delta-potential. 
We also assume that U 5 is low enough so that Hi can be 
treated perturbatively within the linear response theory. 

As in classical physics, the drag force F acting on the 
fluid is linked to the mean dissipated energy per unit time 
by the definition 


(E) = -F-v. (3) 

In a superfluid, displaying frictionless flow, we expect no 
energy dissipation and the drag force vanishes in a certain 
range of velocities. Notice that energy dissipation can be 
directly probed in experiments [2711311116]. Using linear- 
response theory and the fluctuation-dissipation theorem, 
in the thermodynamic limit, the most general form of the 
drag force in ID is 071 (see appendix [A| for details) 

1 r+00 

dq\U{q)WS{q,qv){l - (4) 

where U (g) is the Fourier transform of the barrier poten¬ 
tial U(x) = and 

V / y TV w 



is the dynamic structure factor m, giving the weight 
of the excitation spectrum. (... )o indicates the quan¬ 
tum statistical average with respect to the unperturbed 
Hamiltonian Hq, and Sn{x,t) = n{x,t) — uq represents 


the local fluctuations of the density operator. The dy¬ 
namic structure factor is linked to the Fourier transform 
of the density-density linear response function by the 
fluctuation-dissipation theorem. In Eq. the dynamic 
structure factor is integrated along the line uj = qv in 
the wavenumber-energy plane: if 5'(g, uj) takes arbitrarily 
small values along this line, i.e. no collective excitations 
are possible, then the drag force vanishes and the flow 
is superfluid. We note that this is a sufficient but not 
necessary condition since there are other factors in the 
integrand in Eq. 

In the following, we shall determine the drag force 
acting on the fluid by computing the time-dependent 
density-density correlation functon, evaluate its Eourier 
transform with respect to time and space to get the dy¬ 
namic structure factor (§, and eventually use Eq. Q. 
We expect the drag force to depend on the stirring veloc¬ 
ity r’, the temperature T, the waist w of the barrier and 
its strength Ui^/w^ the interaction between atoms g, and 
the equilibrium linear density no- The strength of boson- 
boson interactions is expressed through the dimensionless 
parameter 


Eint TTig . . 

W 

^kin TIqH 

ratio of the interaction to kinetic energy of the atomic gas 
in the equilibrium Hamiltonian Hq. It can be fine-tuned 
experimentally by Eeshbach or confinement-induced reso¬ 
nances. We will focus on the strongly interacting regime 
7 ^ 1 . In the limit 7 ^ +oc of infinitely interacting 
bosons we use an exact solution [3H]- At arbitrary inter¬ 
action strength we use a low-energy effective Hamiltonian 
given by the linear Luttinger model. 


III. LIMITING CASE 7 ^ -hoc: THE 
TONKS-GIRARDEAU GAS 


A. Exact Bose-Fermi mapping 


The limiting case 7 ^ -hoc is known as the Tonks- 
Girardeau gas [3H]. A well-known peculiarity of one¬ 
dimensional systems is that in many respects, hard-core 
bosons behave like free fermions [48ll5Q] . This phe¬ 
nomenon, related to statistical transmutation, is due to 
the fact that infinite repulsive interactions impose the 
many-body wavefunction to vanish wherever there is a 
contact between two particles, which is the same con¬ 
dition occuring in a Fermi gas owing to antisymme¬ 
try of the wavefunction. Namely, one can write the 
many-body wavefunction of the hard-core Bose gas 7 /;^, 
where the superscript G stands for “Tonks-Girardeau”, 
in terms of the one of a free Fermi gas denoted by 7 /;^, 
as = n(i,j) sign(a;i - Xj)il}^ {xi,... ,xn)- 

The Bose-Fermi mapping has been also demonstrated at 
finite temperature na, namely, the thermal average of 
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an observable O for the bosonic gas can be obtained as 
(O) = (7) 

n 

where A = e~^'^ ^ is the mapping operator ex¬ 

pressing the Jordan-Wigner transformation from bosons 
to fermions, Wn = are the thermal weights and Z is 

the partition function. From Eq. § follows that the par¬ 
ticle density and density-density correlation functions of 
a Tonks-Girardeau gas coincide with the fermionic ones, 
both at zero and finite temperature, since in this case 
[O, ^4] = 0. In the following we shall exploit this prop¬ 
erty to determine the dynamic structure factor of the TG 
gas. 


B. Dynamic structure factor 


We recall first the zero temperature results. Gomput- 
ing fermionic density-density correlations using Wick’s 
theorem yields, after Fourier transform, the well-known 
result for the dynamic structure factor of the Tonks- 
Girardeau gas miEi] 

S^{q,uj) = [‘^+( 9 ) - w] 0 [w - uj-{q)], (8) 


where 0 is the Heaviside distribution and co’+(g) = 
^{q^ + 2qkF), w_(g) = ^\q‘^ - 2qkF\ are the limiting 
dispersion relations bounding the area where particle- 
hole excitations can occur, due to energy conservation. 
Recalling the Bose-Fermi mapping, it is natural to call 
Fermi wavevector the quantity kp = Trno, with no = 
No/L. 

At finite temperature, the linear density-density re¬ 
sponse function for the Tonks-Girardeau gas in the 
Fourier space is given by Lindhard’s expression [5^ 


Xnn(^,^) 


1 Tipik)—npikq) 

L ^ hjjj e{k) — e{k + g) + i0+ ’ 


(9) 


where npik) = is the Fermi-Dirac distribu¬ 

tion and e{k) = the dispersion relation. Using the 
fluctuation-dissipation theorem, we deduce from Eq. © 
the expression of the finite-temperature dynamic struc¬ 
ture factor in the thermodynamic limit: 

S^{q,co)= ( 10 ) 

J —oo ^ ^ 

where ujq{k) = \ [^{k + g) — ^{k)] and 5 is the Dirac dis¬ 
tribution. This is equivalent to 

/ + 00 

dk npik) [1 — npik -h g)] S[(jJ — ujq{k)], (11) 

-oo 

which is easier to interpret physically in terms of particle- 
hole excitations. Either of them can be used to obtain 


S^iq,!^) = 


/ hw-^\ _ ( h:^+e(q) \ 

m "'P K^q/m ) "'P h?qlm ) 

n\q\ 1 - 


( 12 ) 
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FIG. 1. (Color online) Dynamic structure factor Sp{q,uj) of 
the Tonks-Girardeau gas in the thermodynamic limit, in units 
of S^{2kp, 0), for several dimensionless temperatures TjTp — 
0.1, 0.5,1 and 4 in panels a), b), c) and d) respectively, where 
Tp is the Fermi temperature. Frequencies u are expressed 
in units of up = ep/h^ where ep = T?k\/2m is the Fermi 
energy, and wavenumbers g in units of the Fermi wavevector 
kp. Black solid lines correspond to the limiting dispersion 
relations up and oj-, in units of cjf, defining the excitation 
domain at T = 0. 


In order to obtain the finite-temperature dynamic 
structure factor, we have first found the temperature de¬ 
pendence of the chemical potential as detailed in Ap¬ 
pendix The results are shown in Fig. S{q^uj) is 
very sensitive to temperature. Non-vanishing contribu¬ 
tions spread beyond the particle-hole excitation spec¬ 
trum boundaries at finite temperature, since the dynamic 
structure factor includes thermally-activated excitations. 
The latter can even occur at cj < 0, meaning that collec¬ 
tive excitations can be emitted. The quasi-linear shape 
of the spectrum near the origin and the umklapp point 
(2/ci7,0) at T = 0 fades out at temperatures larger than 

O.lTp, where Tp = 2mkB Fermi temperature. In¬ 

creasing the temperature contributes to breaking down 
the imbalance between the cj > 0 and cj < 0 domains, 
see panel d) of Fig. 


C. Drag force 


Using the exact expressions for the dynamic structure 
factor, the next step is to compute the drag force. We 
start recovering the known result in the case re = 0, T = 
0. Combining Eq. @ and Eq. ([T^, we find, in agreement 
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with [37]: 


F^{v) 


2U^nQm 

Ji? 


Q(v — vf) H- Qivp — v) 

vf 


(13) 


where vf = is the Fermi velocity. The drag force 
is linear with the barrier velocity \i v < vf and satu¬ 
rates if V > Vf- As we will see below, this saturation 
is an artifact due to several theoretical simplifications. 
Equation (13) shows that the drag force never vanishes if 


the velocity of the perturbing potential is finite, meaning 
that energy dissipation will occur as long as a barrier is 
driven across the fluid. Hence, the Tonks-Girardeau gas 
is not superfluid according to the drag force criterion. 

We next generalize our calculations to the case of a 
finite laser waist, still assuming that T = 0. The Fourier 
transform of the potential in Eq. reads 


U{q) = Fe- — 


(14) 



FIG. 2. (Color online) Drag force in units of F^{vf) 

in a Tonks-Girardeau gas as a function of the dimensionless 
barrier velocity vjvp- Solid lines stand for a dimensionless 
waist wkF — 0, dashed lines for wkF — 0.5 and thick lines 
for wkF — 1. For a given set of curves, temperature increases 
from 0 to O.ITf to 0.5 Tf from top to bottom, in black, red 
and blue respectively . 


so that the drag force is readily obtained as 




^+oo 

/ 

Jo 


dge 4 qS{q^qv). 


(15) 


IV. FINITE INTERACTION STRENGTH 7: 
THE LUTTINGER LIQUID MODEL 

A. The model 


For the Tonks-Girardeau gas, this yields the analytic ex¬ 
pression at T = 0, 


FSiv) 

FG{vf) 



wkF 




where h{x) = due~'^^. 

We next discuss thermal effects on the drag force. We 
first treat the limit case of an infinitely thin barrier (w 
0) at finite temperature. The drag force reads 


_ 1 [T^ rPrnv^/2 

FG{vf)~ 2]I^ Jo yi(e^-/3^(T) + 1 )• 


(17) 


The integral can easily be evaluated numerically. Notice 
that, at very low temperatures, it is very close to the 
Fermi integral F_i/2{/3eF) [SSI (54]. The most general 
case of finite waist and temperature is obtained by in¬ 
serting Eqs (12) and ( [l4| ) in Eq. @ and then evaluating 
it numerically. Figure shows the drag force as a func¬ 
tion of the potential barrier velocity for several values of 
the temperature and waist. As a main result, thermal ef¬ 
fects broaden the curves around the Fermi velocity, and 
at low velocity the drag force remains linear. If re = 0, 
F{v <C Vf) — ^ thus measuring the slope near 

the origin yields the chemical potential. The finite width 
of the laser beam should be compared to the interparti¬ 
cle distance: the drag force decreases at increasing wirno. 
The saturation of the drag force disappears at finite bar¬ 
rier width. 

Now we proceed to study the case of finite but large 
interaction strengths. 


The Luttinger liquid model was introduced to de¬ 
scribe interacting fermions in ID, where the Fermi liq¬ 
uid paradigm breaks down [HISSID]. Owing to statis¬ 
tical transmutation, the same approach handles as well 
one-dimensional interacting bosons [HISS]. At T = 0, 
Luttinger liquids belong to the universality class of ID 
systems with gapless, linearly-dispersive excitations. The 
Lieb-Liniger model belongs to this class at low enough en¬ 
ergy, since its dynamic structure factor can be linearized 
around the origin and the umklapp point in the {q,oj) 
plane. Therefore, we describe the bosonic fluid by the 
effective Tomonaga-Luttinger Hamiltonian m 




hvs V 
27r Jo 


dx 


K{d,cpf + -{djf 


(18) 


where the superscript LL will hereafter denote quanti¬ 
ties computed for a Luttinger liquid. In Equation (18), 
(j){x) is the phase field in the phase-density representa¬ 
tion of the bosonic field operator '^(x) = 

0{x) is a field related to the number of particles between 
the origin and the position x and whose derivative has 
a peak whenever a particle is encountered, satisfying the 
commutation relation [dx0{x)^ (l){x')] = i7rS{x — x'). The 
isothermal sound velocity Vg and the dimensionless Lut¬ 
tinger parameter K = kmra^jmvs are taken as an input of 
the theory. They can be extracted experimentally since 
they are linked to measurable quantities (compressibility, 
density), or from a microscopic model whose low-energy 
limit is the Luttinger liquid we consider. 

For ID bosons described by the Lieb-Liniger Hamil¬ 
tonian we have obtained the Luttinger parameters 
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FIG. 3. (Color online) Luttinger parameters Vs in units of 
vf (blue circles) and K (dimensionless, purple squares) as a 
function of the dimensionless interaction strength 7 at T = 
0. The solid lines are the asymptotic expansions at large 
interactions. 



FIG. 4. (Color online) Temperature dependence of the Lut¬ 
tinger parameters Vg in units oi vp (blue circles) and K (di¬ 
mensionless, purple squares) in the Tonks-Girardeau regime. 
Solid lines are low-temperature asymptotical expansions, ex¬ 
tracted from the Sommerfeld expansion. 


by numerically solving the Bethe-Ansatz equations (cf. 
Appendix both at zero m and finite temperature 
|5Q| . Results at zero temperature are shown in Fig. 
together with the asymptotic expansion at large interac¬ 
tions, (7) ^ 1 - i + ^ + (^2 _ 6 ) ^ - ( 2^2 -3)1^ 

|63| . The Tonks-Girardeau gas limit 7 ^ +oc treated in 
the previous section corresponds to the values K = 1 and 
Vs = Vf in the Luttinger liquid description at zero tem¬ 
perature. For bosons with repulsive interactions one has 
iF > 1, whereas for repulsive fermions K < 1. The case 
K = 1 corresponds to both infinitely-interacting bosons 
and free fermions, owing to the Bose-Fermi mapping. 

The Luttinger parameters generally depend on temper¬ 
ature. To illustrate this, we compute the temperature- 
dependence of the sound velocity in the Tonks-Girardeau 
regime as extracted from the static structure factor m 
S{q) = dujS{q^uj) according to the compressibility 
sum rule m 




q^O 


mvt 


where S^{q) is evaluated numerically using Eq. (12) 


and the temperature dependence of the chemical po¬ 
tential (see again Appendix H- Then we use the re¬ 
lation Kvs = , stemming from Galilean invari¬ 

ance, hence also valid at finite temperature, to obtain 
the Luttinger parameter K. We also extract Vg ana¬ 
lytically from the Sommerfeld expansion of the chemi¬ 
cal potential at low temperature, yielding for T Tp' 

- l-fj (^) -W (^) ■ Our results are shown 

in Fig. 1^ Note in particular that for T ^ 0.5T;?, one has 
K > 1. Although the effect is small at low tempera¬ 


tures, as we shall see in Sec. |IV D| below it is important 
to include the temperature corrections in the Luttinger 
parameters in order to find agreement with the Tonks- 
Girardeau exact solution. 

We proceed by diagonalizing the Luttinger Hamilto¬ 
nian (p! 8 |) , using a mode expansion of 0 and (j) over bosonic 


fields obtained for periodic boundary conditions m 


9{x) = 9o - 


and 


^(Ar-7Vo) + E 

q^O 


ttK 


2qL 


e-^{e^i^bg + h.c.) 

( 20 ) 


(t>{x) = (t>o + — J +y^sign(g) 


2KqL 


e|g| 


2^""6g + /l.C.), 


( 21 ) 

where h.c. means hermitian conjugated, with hq a 
bosonic field operator satisfying the commutation rela¬ 
tion = ^q,q' and e is a model-dependent high- 

momentum cut-off. The zero-mode terms in Eqs (20) 
and contain respectively the particle number opera¬ 
tor N and the topological number operator J, which we 
shall drop as we do not describe angular-momentum car¬ 
rying states, as well as their conjugate zero-mode fields 
6 >o and 0o- The latter, though may play a role in finite- 
size systems [64l[65], do not contribute to the correlation 
function of interest here. Once diagonalized, the Lut¬ 
tinger Hamiltonian (H 8 |) reads 


H 


LL 


= E ^i9)b\hq 

q^O 


( 22 ) 


with uj{q) = \q\vs^ i.e. the bosonic field obeys a linear 
dispersion relation. These phonons are collective oscilla¬ 
tions of the phase and density fields and well describe the 
excitations of the bosonic fluid at low momentum. We 
proceed to compute the density-density correlations, the 
dynamic structure factor and the drag force. 


B. Density-density correlations 

We start by recovering the zero temperature result. 
Using the bosonized expression of the density opera¬ 
tor: n{x) = ^dx0{x) X]m=-oo [61] and 
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the mode expansion ( 20 ) we obtain in agreement with 
iniss], up to leading orders in the index m, 


{6n{x^ t)Sn{0,0 ))^^ 


K 


Ak‘^ {x — Vgt + ieY {x + Vgt — ie)^ 
cos(2/ci7x) 




riQ [{x — Vgt + ie){x + Vgt — ie] 


iK 


, (23) 


where Ai is a non-universal parameter. We can express 
it as Ai{K) = 2(eno)^^ where the cut-off e depends on 
the interaction strength 7 and thereby on the Luttinger 
parameter K. As we shall see below, the first term is 
related to the dynamic structure factor near the origin 
in the (g, 00 ) plane, while the second one has a finite con¬ 
tribution in the vicinity of the umklapp point ( 2 /ci?, 0 ). 
The algebraic decay of the correlations is a signature of 
quasi-long-range order. We have kept only the two lead¬ 
ing terms since larger orders decay with an increasingly 
large coefficient Km? ^ with m a non-negative integer. 
This expression does not depend on the statistics, except 
through the values of K and which span different 
ranges for repulsive fermions and bosons. 

Using the same approach, we generalize the previous 
result to finite temperature (details are found in Ap¬ 
pendix 


{6n{x^ t)Sn{0, 0))t^ 


K TT^ 


sinh^ 


7r{x — Vst) 
L'X' 


T + 



1 

sinh^ 

7T(x-\-Vst) 

L'jp 


{^sinh 


2 cos(2/ci?x) 


7T(x — Vst) 

Lt 


sinh 


Tv{x + Vst) 1 I 



where Lt = Phvg is a thermal length. This expression 
is valid in the limit x ^Vst^Vgt^ Lt ^ e, and agrees with 
[62] . Here again we have kept only the two leading terms 
in the dynamic factor because of the exponential decay 
of higher-order terms with a Km^ exponent. 


C. Dynamic structure factor 


At T = 0, the Fourier transform of Eq. ( |23[ ) yields the 
dominant terms of the dynamic structure factor. Since 
the latter is symmetric with respect to the cj-axis, we 
write the result for g > 0 


S^^iq,Lo) = K\q\5[uj - uj{q)] 

+Bi{K) [a;2 - {q - 2kFfv^sf~' " k " 

= 57 {q, Lo) + S^^{q,u), (25) 


where Bi{K) = ^ 2 notl)nLi} is a non-universal 

coefficient. In the dynamic structure factor, Sq displays 



FIG. 5. (Color online) Definition domain of the dynamic 
structure factor at T = 0 in the plane (g,ct;) in units of 
{kF^LUp). We superimposed the result for a Tonks-Girardeau 
gas to the result in the Luttinger liquid framework for di¬ 
mensionless parameters K = 1, Vg/vp — In the latter, 
the area consists in a line starting from the origin, and the 
area included in the triangle starting from the umklapp point 
(0,2/ci?). The upper energy limit of potential validity of the 
Luttinger liquid model is approximately given by the dashed 
line. 


a sharp feature in correspondence to the linear dispersion 
uj{q) = qVg. There are also two linear limiting dispersion 
relations described by Si , symmetric with respect to the 
q = 2kF line, forming a triangular shape at the umklapp 
point (2/ci7,0). The slopes of the limiting dispersions in 
So and Si depend on the interaction strength via the 
interaction-dependent sound velocity. Hence, measuring 
them for the Lieb-Liniger gas is a way to determine Vg. 

In Fig. we plot the definition domain of the dynamic 
structure factor of a Luttinger liquid with parameters 
K = 1 and Vg = vp, superimposed on the result found in 
the previous section for a Tonks-Girardeau gas The 
comparison shows the domain in the (q^uj) plane where 
the Luttinger liquid description is valid: namely, at low 
energy, at small g and for g 2k close to the umk¬ 
lapp point, where the curvature of the dispersions uj±{q) 
can be neglected. By comparing Eq. ( [^ with Eq. (H, 
we notice that S^^{q^ uj) close to the umklapp point does 
not reproduce the 1/g behavior found in the exact model, 
yet there is quantitative agreement with less than 10% 
error all the same, provided that uj < O.SScj^. This re¬ 
sult coincides with the limiting domain found graphically 
on Eig. I^by looking at the curvature. In this analysis, 
the cut-off e is estimated by stating that the matching is 
optimal on the g = 2kF line, yielding Bi{K = 1) = 
and e{K = 1) = respectively. 

We derive next the dynamic structure factor of a Lut¬ 
tinger liquid at finite temperature. The term SQr^{q,uj) 
is more easily computed in one step rather than by evalu¬ 
ating the real-space density-density correlation functions 
as an intermediate result. Bosonizing from scratch at 
finite temperature, we find the first contribution to the 
dynamic structure factor (details can be found in Ap- 
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pendix [E| ) 

- "(9)1 + + "(9)|}-(26) 

By comparing with the exact results found in Sec. |mB] 
we can now discuss the issue of the temperature range 
over which a Luttinger liquid may yield results close to 
those found with the Lieb-Liniger model. Equation (26) 


shows that in the Luttinger liquid framework the dy¬ 
namic structure factor at finite temperature remains lin¬ 
ear near the origin, and has two branches. The one 
corresponding to u; < 0 has a lower weight than the 
cj > 0 one and disappears at T = 0. This behavior is 
observed in the Tonks-Girardeau limit at low tempera¬ 
ture (see again Fig. [^, yet the linear Luttinger liquid 
theory does not predict the thermal broadening of the 
dispersion relation, which becomes more and more rele¬ 
vant at increasing temperature, nor the broadening and 
the curvature due to non-linearities. We estimate that 
in the infinitely interacting regime, the Luttinger liquid 
predictions for the dynamic structure factor around the 
origin are satisfactory for temperatures lower than ap¬ 
proximately O.lSTi?. At larger temperatures, it is not 
relevant to linearize the Tonks-Girardeau dynamic struc¬ 
ture factor around the origin because the broadening is 
too pronounced. Furthermore, combining Eq. (26) with 


the linearization at small q of the exact dynamic struc¬ 
ture factor, we verify that the temperature dependence 
of the sound velocity is small at low temperature: for 
T < 0.15Tp, we find 0.97 < Vg/vp < 1.00, in agreement 
with Fig. 

The expression of the backscattering term 
in the dynamic structure factor is obtained from the 
Fourier transform of the density-density correlations. Us¬ 
ing the property [67] r(— -h^) = (2i)^+^^r(l -1- 

sin^(t^), > 0,i?e(x — yzi) > 0, 

where T is the complex Euler Gamma function, after 
some algebra, we obtain (for more details, please refer 
to Appendix [E|) 

Siril > 0,w) = C{no,Vs,e,T)e^ 


B 


B 


K jn, . , K ./3n, . ; 

y+y “ 

K pn, ^ , K ph, , ; 

y +~ y “ ~ 


, (27) 


where C{no,Vs,e,T) = q = q - 2kF 

and B[x,y] = is the Euler Beta function. The 

case K = 1 can be computed separately, leading to 

Sit{q > 0,w)|if=i 


2u,, 


cosh 


^{oj + qvs) 


cosh 




in agreement with the general case since r[z] = r[z] and 

|r[| + *i/]r = 


cosh[7ry] ' 



2 q/kp ^ ^ ^ q/kp ^ 


FIG. 6. (Color online) Dynamical structure factor in units 
of 5'(2/cf, 0)t=o at T = O.lTp in the plane {q,uj) in units of 
(kF,(^F) in the vicinity of the umklapp point, as predicted 
for a Tonks-Girardeau gas (left panel) and a Luttinger liq¬ 
uid for dimensionless parameters K = 1.005, r’s/r’F = 0.995 
(right panel). The exact temperature-dependence is quite well 
reproduced by the Luttinger liquid model, differences come 
mostly from the non-linearities which are not taken into ac¬ 
count in the Luttinger liquid framework. 


As we did for Sq, for the backscattering term Si too we 
can assess the regime of validity of the Luttinger liquid 
description at finite temperature by comparing with the 
Tonks-Girardeau results. Fig. |IV C| shows the dynamic 
structure factor of a Tonks-Girardeau gas at tempera¬ 
ture T = O.lTp and the one computed for a Luttinger 
liquid with the appropriate Luttinger parameters at this 
temperature. The cut-off e(O.lTp) is chosen so that the 
matching is optimal on the q = 2kF line, i.e. they coin¬ 
cide exactly at low energy. Around the umklapp point 
we find that the Luttinger liquid reproduces quite well 
the exact thermal broadening, provided that the energy 
is low enough so that non-linear effects can be neglected, 
as was already the case at T = 0. 


D. Drag force 


Once the dynamic structure factor is known, we can 
compute the drag force. First, we address the limit T = 
0, re = 0, which yields 


tj 2 n+oo 


dqqS^^{q,qv) 


C/2 B^{K) yPT{K){2kFVsf^ 
2TTh vl T{K + 1/2) 


/ / s,2if-l 

{v/Vs) 

l-{v/Vaf 


(29) 


in agreement with the result found in [28| in the limit 
v/vs 1: at low velocities, the drag force scales as 
a power law A comparison with the Tonks- 

Girardeau result at iC = 1 leads to the determination 
of the constant Bi{K = 1) = 

Then, we generalize the expression of the drag force 
to finite laser waist w. For the case K = 1 we find the 































FIG. 7. (Color online) Drag force in units of as a 

function of the velocity v (in units of vf), as predicted for 
a Tonks-Girardeau gas (dashed lines) and a Luttinger liquid 
model at dimensionless parameter K — 1 (solid lines), at 
T = 0. Thin black curves correspond to a dimensionless waist 
le/cF = 0 and thick red curves to a finite waist wkF = 0.5. 


analytical expression 


:^LL 




U'jnom 

2tfik‘pw‘^ 




(30) 


We use Eqs (13),(16),(29) and (30) to plot the curves in 
Fig.0 showing that the Luttinger liquid model is able to 
reproduce exact results in the Tonks-Girardeau regime 
at low velocities. The drag force is all the better approx¬ 
imated as the potential is wide. It is always linear near 
the origin, with a slope depending on w. 

At arbitrary interactions, the expression of the drag 
force for the case of a finite-width potential is given by 


,LL. B^{K) ^T{K){2kpv,f^ {v/v,) 


2K-1 


Fi;\v) = 


+00 


27rh V{K F 1/2) 

2/c+l 


(—1)^ / wkp 


wkp^ k\ 

^ k=0 


.1 + -. 


Fi -l-2k,K;2K;- 


1 - (v/Vs) 
2v 


Vs-V 


1 K 


(31) 


where 2 F 1 is the hypergeometric function. We have ver¬ 
ified that for wkp % 1, truncating the sum at low orders 
yields a very good accuracy. 

The effect of temperature on the drag force is ob¬ 
tained by integrating numerically Eq. with the input 
of Eqs (26) and ([^. We have plotted in Eig.j^the drag 


force at T = 0 and finite temperature as a function of the 
velocity for a Tonks-Girardeau gas as obtained from the 
exact solution and the Luttinger liquid approach, where 
the value of the cut-off e is chosen by enforcing that the 
results should coincide at the origin. We note that we 
find an excellent agreement at lox velocities by taking 
into account finite-temperature corrections of the LL pa¬ 
rameters while the use of zero-temperature values would 
yield a less precise agreement (not shown) The predic¬ 
tions of the Luttinger liquid and the Tonks-Girardeau 
gas start differing at the same velocity whether tempera¬ 
ture is taken into account or not, essentially due to non- 
linearities of the Tonks-Girardeau dispersion. 



FIG. 8. (Golor online) Drag force in units of F^{vf) as a 
function of the velocity v (in units of vf), as predicted for 
a Tonks-Girardeau gas (dashed lines) and a Luttinger liquid 
(solid lines), at ic = 0. Thin black curves correspond to T = 0 
and thick blue curves to T = O.lTp. 


As a main result, we have shown that in the regime of 
very large interactions, the Luttinger liquid theory is able 
to reproduce the exact results of the Tonks-Girardeau gas 
in terms of dynamic structure factor around the umklapp 
point and drag force at low velocities, even for a potential 
barrier with a finite width. This allows us to use the 
Luttinger liquid theory to predict the generic behavior 
of the drag force at large to intermediate interactions, 
thus complementing the Bogoliubov approach at weak 
interactions. 


V. SUMMARY AND OUTLOOK 

We have used the concepts of dynamic structure factor 
and drag force to explore theoretically the superfluidity 
of a system of strongly-interacting bosons in ID, stirred 
by a Gaussian laser beam. We have studied the limit¬ 
ing case of infinite interaction strength using the exact 
Bose-Fermi mapping, and compared it to the predictions 
of the linear Luttinger liquid model. We have obtained 
various analytical expressions generalizing known results 
to finite temperature and finite laser beam width. The 
Luttinger liquid model predictions are limited by the non¬ 
linearity of the real physical system which are not taken 
into account within this model. As seen in the exact 
Tonks-Girardeau solution, at higher energies or at in¬ 
termediate wavevectors q kp^ beyond-linear Luttinger 
liquid effects will appear [SEHZD- These have been ex¬ 
perimentally observed [32]. Yet, our work shows that the 
effects of temperature on the dynamic structure factor 
around the umklapp point {q = 2kF^uj = 0) and on the 
drag force are well taken into account. Our results lead 
to precise estimates for the low-energy behavior of the 
dynamic structure factor at finite temperature. General¬ 
izations of the drag force to a finite-size potential barrier 
show that this parameter has a dramatic impact, which 
is well reproduced by use of the Luttinger-liquid model 
to treat the Tonks-Girardeau regime. We conclude that, 
provided that the temperature is low enough and the ve¬ 
locities small, the Luttinger liquid theory can be used 
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to test the superfluidity of the Lieb-Liniger gas accord¬ 
ing to the drag force criterion. Moreover, our results are 
expected to describe a wide range of systems whose low- 
energy description is a Luttinger Liquid, for instance, the 
Calogero-Sutherland model mizg. The barrier poten¬ 
tial we have chosen may be realistic to describe future 
experiments, yet to be even more realistic, finite size ef¬ 
fects e.g. in ring traps, should be taken into account, as 
well as the possible local non-homogeneity of the fluid, 
e.g. in harmonic traps m- 
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Appendix A: Drag force in the linear response 
theory 

In quantum physics, if a system is described by a 
Hamiltonian H = Hq ^ Hp^rt = + / d^rA{r)U{r,t), 

where A is a linear operator and U a weak local pertur¬ 
bation, then the average of any observable B coupled to 
A in the presence of the perturbing potential reads 

/ + 00 p 

dt' d'^rU (/, t')xBA{r-r', t-t') (Al) 

-oo j 

to first order, defining the B — A linear response function, 
denoted by Xba- Using the interaction representation 
(denoted by a subscript I) and Liouville equation for the 
density matrix, one can show that if the system without 
perturbation is Galilean-invariant, then 



FIG. 9. (Color online) Chemical potential fi of the Tonks- 
Girardeau gas in units of the Fermi energy cf as a function 
of the reduced temperature T/Tp, where Tf — f\^k^ j{2mkB) 
and kp = ttuq. The red solid line stands for Sommerfeld’s low- 
temperature expansion to order 4: ji/ep — 1 -\- /TfY + 

^^(T/Tp)^, blue points are the numerical result. The posi¬ 
tion (To/Tp,0) of the isolated point was computed analyti¬ 
cally, stating that ji should vanish at this temperature, yield¬ 
ing To/Tp = 4 /7r[(v^ — l)^(l/2)]^ ~ 3.48, where ( is the 
Riemann zeta function. 


parts: Xnn{q,A = x'nniq^A + The real part 

is an even function of its arguments while the imaginary 
part is odd. Together with the fluctuation-dissipation 
theorem, this yields: 






(A5) 


One then specializes to the ID case and uses the fact that 
if the perturbing potential is a function of x — vt, then 
its Fourier transform contains S(cj — qv) due to energy 
conservation, to find, eventually, Eq. with U{q^uj) = 
2ttU {q)6{uj — qv). 


XBA{r-r' = {[Bi{r-r',t-t'),A])Q (A2) 

Let A = B = n the density operator. The ensemble 
average dissipated energy per unit time reads: 

{E) = j (VrU{r,t)n{r,t)) ■ (A3) 

Using the continuity equation and the definition 
of the density-density response function, with the 
convention for the Eourier transform f{q,oj) = 
f {x ^ t), we obtain: 




(27r) 


{E) = J d'^r J J 

We then compute the time average, use the property 
U{—q, —uj) = lB^{q^uj) and decompose the Eourier trans¬ 
form of the response function into its real and imaginary 


Appendix B: Chemical potential of the 
Tonks-Girardeau gas at finite temperature 

In Eig. we show the temperature dependence of the 
chemical potential of the Tonks-Girardeau gas, which co¬ 
incides with the one of a one-dimensional ideal Eermi 
gas. We note in particular that at low temperature, due 
to the reduced dimensionality, the chemical potential in¬ 
creases with temperature, at difference from the three- 
dimensional case. 


Appendix C: Solving the Bethe Ansatz equations of 
the Lieb-Liniger model 

In [50] the Lieb-Liniger Hamiltonian for N bosons is 
written as 


N 


i=l 


i>j 


(Cl) 
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FIG. 10. (Color online) Chemical potential /i in units of cf 
as a function of the inverse dimensionless interaction strength 
1 / 7 , at T = 0 (blue circles) and T = O.ITf (red squares). 


The Bethe Ansatz solution for the equation of state at fi¬ 
nite temperature in the thermodynamic limit is obtained 
by solving a system of three coupled equations: 


e{k) = —fi^c^T) + 


Tc /'+°° dq 
TT y_oo c2 + (fc - q)2 


In 


11 + exp 



27Tp{k) ^ 1 + exp 
and 


e{k) 




— 00 


(k — q)‘^ 


, (C3) 


/ + 00 
-00 


dk p{k) = no, 


(C4) 


where the k are quasi-momenta, p is the quasi¬ 
momentum distribution function and exp[e(/c)/T] = pf^/p 
with ph the quasi-hole distribution function. With the 
correspondence c = yno, we extract p. (j^T/Tf) /ep from 
this set of equations. To find the chemical potential at a 
given reduced temperature and interaction strength, we 
proceed as follows. We g uess the value of p and self- 
consistently solve Eq. (C2) by iteration to find the func¬ 


tion e{q). On ce e is found with enough accuracy, we inject 
it in Eq. ( |C3[ ) and find p by iteration. When fluctuat ions 
become low enough to be neglected, we check if Eq. (C4) 


is verified up to the chosen accuracy. Most probably, the 
initial guess for p is not precise enough, so this procedure 
is repeated as many times as needed within a simplex al¬ 
gorithm to find p with the chosen accuracy. A few steps 
are enough to find p/ep with less than 0.5% error. Our 
results are shown on Eig. (10) and Eig. 0- The chem¬ 
ical potential depends both on the temperature and the 
interaction strength. When the latter is decreased, the 
chemical potential increases, while it increases with the 
reduced temperature. 



FIG. 11. (Color online) Chemical potential p in units of cf 
as a function of the dimensionless temperature T/Tp for a 
Tonks-Girardeau gas (blue circles), for a dimensionless inter¬ 
action strength 7 = 1000 (red squares) and 7 = 100 (yellow 
diamonds) respectively. 


Appendix D: Details of the calculation of the 
density-density correlations in the Luttinger liquid 
framework at finite temperature 


Here we sketch the derivation of Eq. (24). Using 
n{x) = -dx0{x) ^2im[9(x)-\-kFx] havc 


{n{x,t)n{0,0)) = -^{d:c0{x,t)dxd{0,0)) 

1 1 

_j \ ^ ^2irnkFX ^^2irn9(x,t) ^2irn'9(0,0) (El) 

7r2 

m,m' = — l,j^(0,0) 


which will lead to the zero-order and first-order terms 
respectively. To find the first-order term we in¬ 
troduce a generating function: Gm,m '= 
pimeix,t)pim'e(0fi)_ To evaluate {Gm,m>{x,t;0,0)) we use 

the identity valid for any pair of 

operators A and B commuting with their commutator, 
and the property (e^) = \ valid for any linear op¬ 

erator A, together with the mode expansion of the field 0 
(20) and the bosonic commutation relations. After some 
algebra, we find 


g-I I -1}] 

g- Eq^o I ^ I [(m+m')‘^+2mm' {cos[qx-uj(q)t]-l}]nB (q) ^(D2) 


where nB(^) = is the Bose-Einstein distribution 

for the phonons. In the thermodynamic limit the non¬ 
vanishing contributions are those where m = —m!. We 
find 


(^5 ^5 O 5 0 )) — 

^—2Km? cos(gfc)+2[l—cos(gfc) cos(qvst)]n b (q)} 

= (D3) 

To compute F(x,t), one can rewrite ^ = 



















11 


f 

Jo 


Y —>-+oo 


dye yielding: 


^2F(x,t) _ - ie){x - Vgt + ie) 


+00 r 

n 

n=l 


1 + 


The property 


(a; + Vst)'^ 
(e + uLt)^ 

2 


rw 


r{x-iy) 


+00 r 

n 

m=l 

r + OO 


1 + 


(x - Vsty 


(e + mLTY 


(D4) 


- Tf^ 

“ llfc =0 


1 + 


(a:+/e )2 


^2F{x,t) = 


yields: 

1 


_l_ {x+Vsty 2 _j_ 




L'jp 




(D5) 


then in the limit Vgt^x ± Vgt^ Lt ^ e, the properties 
r{a;) ^^^0 I and \T{iy)\ = yield Eq. (|^ after 

a few rearrangements. 


Appendix E: Details of the calculation of the 
dynamic structure factor in the Luttinger liquid 
framework at finite temperature 


yielding Eq. (26). To prove Eq. (27), we start from the 
second contribution in Eq. (24) and compute its Eourier 
transform, yielding 


(x/i(a)/i( 6 ), (E5) 

where Ii{x) = sinh(i 4 )“^, a = [^ + 

{q - 2kFVs)] and b = ^[uj - {q - 2kFVs)]. Using 
d^e-^^^sinh(^)-^ = 2^-ir(^4^)r(i-A)/r(i + 
^^ 2 ^ ) treating the branchcut carefully, we find 


^IT {Q > O5 <^) — 


phvs 


2{1-K) 2 

2Va 


r( 


K -|-i(2 ^ 
2 ) 


—iKiv 


p _ K—ia 

r( 


2 

^ K-ih 


r(^) 


) 


r(i- 


K+ih 

2 


+ e 


iKn 


r( 


K±ib^ 


r(i- 


K-ib\ 


,(E 6 ) 


then after some algebra, using twice the property 
r( 2 ;)r(l - z) = we finally obtain Eq. (27). In 

the case A = 1 , the property dx = 7 r[i + 

cotan (TT/i) ], 0 < Re{fi) < 1 yields Eq. (28) in a few lines 
of algebra. 


Here we sketch a derivation of Eqs. (26) and (27). To 
prove Eq. (26), we split 5'^^(g,ci;) into two contibutions: 


S^^{q,uj) = Sy'^-\q,Lo) + Sy^{q,uj) (El) 


Appendix F: Details of the calculation of the drag 
force in the Luttinger liquid framework at finite 
barrier width 


where the first part is the result at T = 0 and the second 
is a purely thermal part. We focus on 5'o cc), which 

is more easily computed starting from an intermediate re¬ 
sult in the calculation of the density-density correlations. 
In order to compute 


oLL,Te X K 
So {q,uj) = ^ 


/ Yoo pYoo 
00 J —00 


[ dg|g|nB(g) (E2) 

Jq^O ^ ^ 


we perform the change of variables u = x — Vgt and v = 
X yvgt. After some algebra we find 


^0 "'’^(5, - a;(g)] +(5[a; + Y<l)] XE3) 

thus 


SoriQ^Y = K\q\5[uj - uj{q)] 

g/3Rw(!^L 1 ~ I (E4) 


Here we derive Eq. (jsT]). We need to evaluate: 
/y dg q{q - q-Y~^{q+ - qY~^ exp We split 

it into two parts using: q = g —g_-hg-, then expand the 
exponential as a power series: exp(x) = ^ and 

expand once again according to g^^ = (g —g_-hg-)^^ = 
Ylm=o (m)(^ “ q-)'^q^_^~'^. To eliminate the remaining 
integrals, we use the property dx{x—a)^~^{b—xY~^ = 
( 6 —z/), where B is the Euler Beta function, 
then after a resummation and using Pascal’s triangle, we 

are left to evaluate W = Ylm=o ( i-i]vs ) B{K+m, K), 
which can be interpreted as the sum of a series with all 
its terms equal to 0 from rank 2k d-2 on. We express it in 
terms of the hypergeometric series 2 E 1 , which converges 
since— 2 /c — l is a negative integer: 


1 T = 


r(A-hi/ 2 ) 


-l-2k,K: 2K: 



This readily yields Eq. (31). 
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